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ABSTRACT 

q . The Westerbork Synthesis Radio Telescope finds a weak 21cm line emission 

feature at the coordinates (RA-Dec-velocity) of the sub-Damped Lyman-a ab- 
sorber observed at z a b s = 0.00632 in the spectrum of PG1216+069. The emission 
feature, WSRT-J121921+0639, lies within 30" of the quasar sightline, is detected 
^ at 99.8% (3a) confidence level, has Mhi between 5 and 15xlO 6 M , and has ve- 

locity spread between 20 and 60 km s _1 . Other HI emitters in the field include 
VCC297 at a projected distance of 86/if 5 1 kpc and a previously unreported HI 
cloud , WSRT-J121919+0624 at 112/if B 1 kpc with M m ~ 3xlO 8 M . The opti- 
cally identified, foreground galaxy that is closest to the quasar sightline appears 
to be VCC339 (~L*/25) at 29/i 75 1 kpc and velocity offset 292 km s" 1 . A low 
surface brightness galaxy with the HI mass of the sub-DLA absorber WSRT- 
J121921+0639 would likely have m B ~ 17, and its diffuse optical emission would 
need to compete with the light of both the background QSO and a brighter 
foreground star ~10" from the QSO sight line. 

Subject headings: cosmology: observations — quasars: absorption lines — galax- 
ies: formation - - galaxies: intergalactic medium - - radio lines: galaxies - 
quasars: individual(PG1216+069) 
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1. Introduction 

The goal of clarifying the nature of QSO absorption line systems drives the interest 
in extreme low redshift absorption systems such as the "sub-Damped Lyman-a" absorber 
in the sight line to the z = 0.33 AGN PG1216+069. At a heliocentric redshift velocity of 
1895 km s -1 , the absorber falls in the outskirts of the Virgo Cluster, where its association 
with a particular galaxy host or group/cluster environment should become clear through 
observation. Tripp et al.s' recent, comprehensive study of PG1216+069 (Tripp et al. 2005) 
makes it clear that (1) the metallicity is low, [0/H]f«-1.6, but at a characteristic level for 
high redshift absorbers, (2) there is no detected optical emission from the absorber itself 
nor an obvious association with a specific galaxy in the field, and (3) the HI column density 
through the absorber Nhi ~ 10 19,3 cm -2 is a factor 10 below the Damped Lyman-« (DLA) 
class. While the lower limit to DLA column density 10 2a3 cm -2 was originally set as a 
completeness limit in the early surveys (Wolfe et al. 1986), this column density also takes on 
physical significance in disk galaxies, since it is typical of the threshold column density for 
the onset of star formation (Kennicutt 1989). Zwaan et al. (2005) used an extensive data 
set of 21cm line observations to demonstrate the consistence between QSO absorption line 
statistics and the cross sections that z ~ galaxies present for absorption. Zwaan et al. 
further demonstrate that the decline in metal abundance for low column density absorbers 
agrees with abundance gradients in the outskirts of nearby galaxies. 

QSO absorption statistics indicate through the f(N) distribution (Tytler 1987; Petitjean 
et al. 1993) that populations of low column density absorbers present larger cross sections 
than the DLA class, implying that clouds with Nhi ~ 10 19 cm -2 may be more common 
than the higher Nhi ones in spiral galaxy disks. Radio telescopes routinely observe Nhi 
well below this level, thereby providing the incentive for a WSRT reconnaissance into the 
PG1216+069 field. Kanekar & Chengalur (2005) observed this field with GMRT within the 
same time frame without detecting the new features we report here. 

2. The WSRT Observations 

The Westerbork Synthesis Radio Telescope (WSRT) observed the field centered on 
PG1216+069 in "filler time" scheduled through the interim proposal system. Rather than 
a single long integration in one 12 hour integration, the observation consists of two shorter 
runs. On 2004 August 14, the first session, covering extreme negative hour angles, produced 
4.5 hours on source in the "2 x 144m" configuration (with RT9 to movable telescope dis- 
tances of 144, 288, 1368, and 1440 meters). On 2004 August 27, the second session covered 
4 hours of extreme positive hour angles in the "Maxi-short" configuration (RT9 to movable 
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telescopes distances of 36, 90 133, 1404 meters). The DZB spectrometer analyzed both 
linear polarizations with a 10 MHz (~2000 km s _1 ) band centered on heliocentric redshift 
velocity 1900 km s~ x with 256 spectral channels per polarization to obtain channel spacing 
of 8.3 km s _1 or velocity resolution of 16.5 km s _1 after Hanning smoothing. 

For an east- west array such as WSRT, observing sources close to the celestial equator has 
two important consequences: The negative one is that the telescope resolution in the north- 
south direction deteriorates at low declination 5 as esc 5, so that the best resolution possible 
for PG1216+069 is a north-south fan beam of 108" x 13". The positive consequence is that 
all projected baselines foreshorten at large hour angle h as cos h, so that the apparent lack of 
short spacings in the "2 x 144m" configuration data does not damage the surface brightness 
sensitivity to extended emission. Short spacing sensitivity is limited to >25m projected 
spacings in order to avoid shadowing of one dish by its neighbor, causing the observation to 
be insensitive to smooth structure on angular scales larger than ~15 arcminutes (115 kpc). 

Data processing steps included editing, calibration and application of different weight- 
ing schemes in the synthesis imaging software packages AIPS and MIRIAD. This produced 
spectral line data cubes with a range of resolution and noise level, summarized in Table 1. 
For each of three weighting schemes, the table lists the angular resolution and noise level as 
measured from the processed data cubes. Assuming that the emission fills the synthesized 
beam solid angle Qb leads to a surface brightness sensitivity (AT# = o"sA 2 /2kB^£) and col- 
umn density sensitivity (ANhi = 1.8xl0 18 AVk m / s ATB cm~ 2 ). The table includes estimates 
of the flux density that a beam filled with Nhi ~ 10 19,3 cm -2 would produce, along with 
the HI mass contained within the beam (Mhi = 2.3xl0 5 AVk m /s5'jy^MpcM©), assuming the 
distance of 26.6 Mpc, which Tripp et al. (2005) adopted. 

Analysis of the continuum image formed from the average of the channels free of HI 
signal found that PG1216+069 had a continuum flux density of 2.7±0.2 mJy at 1416 MHz 
(consistent with the 2.8±0.3 mJy measurement of Kanekar & Chengalur (2005)). The NVSS 
survey (Condon et al. 1998) indicates an upper limit ~1.9 mJy (~2(j), while VLBA obser- 
vations indicated compact structure of 5.7±0.3 mJy in the period 1999-2000 (Ulvestad et al. 
2005). The PG1216+069 AGN appears to be a variable radio source. 

Figures 1 and 2 summarize the observational results. Fig 1 presents contours of the inte- 
gral line emission, obtained from the zeroth moment of the spectral channels after automated 
masking using a smoothed version of the cube. The image is 1 degree on a side, correspond- 
ing to 460 kpc. The half power beam width of the WSRT primary beam is ~35 arcmin, so 
the outer regions of HI map include the first null of the pattern. Table 2 lists the galaxies 
with redshift less than 3000 km s -1 within 25 arcmin of the PG1216+069 sight line; many 
of these are marked in Fig 1. Figure 2 presents spectra for selected areas of the image. 
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3. Discussion 

Figure 1 shows three significant detections for systems VCC297, VCC415, and a pre- 
viously unreported HI emitter, here named WSRT-J121918+0624, centered at coordinates 
12H9 m 17 s .6, 06°24'06" (J2000), close in the sky to VCC327. The HIPASS online database 1 
shows this feature with 5a significance. The HI mass in WSRT-J121918+0624 amounts 
to ~3xlO 8 M contained in a gaussian profile of FWHM ~38 km s" 1 centered on V hd = 
1963 km s -1 , but the 800 km s _1 difference in redshift from VCC327 suggest that this is but 
chance alignment between two unassociated objects. However, the NED 2 indicates that there 
are at least seven objects (NGC4223, VCC287, VCC297, VCC332, WSRT-J121918+0624, 
VCC339, VCC346) with velocities in the 400 km s" 1 range 1870 to 2270 within a 300 kpc 
square area projected on the sky. There is weak evidence that WSRT-J121918+0624 is ex- 
tended toward the NE, in which case the morphology would be typical of tidal debris from 
interactions within the overdensity, similar to the structures seen in Ursa Major (Verheijen 
& Sancisi 2001) and elsewhere in Virgo (Chengalur et al. 1995; Oosterloo & van Gorkom 
2005). 

The catalogued galaxy closest to the PG1216 sight line is VCC329 with a velocity offset 
of 292 km s" 1 relative to the sub-DLA, but there is no hint of an HI signal in the WSRT 
data (see Fig 2). 

Inspection of the spectrum along the line of sight to PG1216+069 reveals a weak emis- 
sion feature at Vh e i = 1905, overlapping the velocity of the sub-DLA absorber. The peak 
strength is 1.7±0.55mJy; for well behaved noise, there is only a 0.14% probability of a 3a 
positive deflection occurring in a pre-chosen pixel by chance. The velocity width is not well 
constrained, and considering widths in the range 20-60 km s -1 would indicate HI masses of 
5 to 15xlO 6 M . The strength of the feature is between 1.5 and 1.7 mJy per beam for all 
three synthesized beams in Table 1, implying that a large fraction of the emission is likely 
to be contained within the 112" x 16" beam. Unless the HI spatial distribution is similar in 
shape to the elongated beam, the HI column density must vary in places to values above that 
specified by the damping wing profile of Tripp et al. (2005). The signal strength measured 
by WSRT indicating 5xlO 6 M for a 16 km s" 1 velocity width is compatible with the GMRT 
upper limit of llxlO 6 M for a 20 km s -1 velocity width and 39.4" x 37.7" angular resolution. 

A rough rule of thumb for estimation of the optical brightness of a late type, gas rich 
galaxy is M B = -20.3 - 2.81og(M H //lO lo M ) (Briggs 1990). An M HI « 1O 7 M translates 
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to M B ~-14.9, with large scatter. With a distance modulus equal 32.1, the integral optical 
emission would be magnitude 17 and possibly much dimmer as would be typical of a low mass, 
diffuse LSB system. The optical emission must compete with the QSO light (magnitude 
15.84, 15.65, and 15.61 in the B, g and i bands, respectively) and a star 10" away that is 
tabulated at g = 14.9 (according to Kirhakos et al. (1994) derived from an image in which 
it is saturated). The APMCat 3 evaluates the sum of the star plus QSO as R = 13.7 and 
B = 14.9. 



4. Conclusions 

The HI emission properties of the PG1216+069 sub-DLA absorber are consistent with 
expectations based on the low column density, low metallicity and absence of optical emission 
from a clear host galaxy. Low column density below the star formation threshold means that 
the lack of star formation locally along the line of sight is no surprise. A low metallicity 
cloud with no connection to a optically identified host galaxy is compatible with the sub-DLA 
originating in a tiny, extreme low surface brightness dwarf of primitive composition or being 
a diffuse cloud in the outer halo of a galaxy such as VCC339 or VCC297. The additional 
presence in this field of the more massive cloud WSRT-J121918+0624 may argue for an 
origin for both features as tidal debris in this overdense volume. (The tranverse dimension 
of the field in Fig. 1 is 460 kpc, thus fitting easily within Local Group, for example.) 

Further progress might be made by (1) longer integrations on this field with the goal 
of tracing low column density structures that may form a more obvious link to the visible 
galaxies in this overdensity and (2) more extensive, deep HI mapping in several other fields 
to identify examples of objects (dwarfs or halo clouds) in more favorable situations (i.e., 
without optically bright QSOs cluttering the field). Such populations might account for the 
sub-DLA and Lyman-limit absorption cross sections. 
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Rene Vermeulen and Tony Foley for coordinating the observations through the Interim Pro- 
posal program. The project was supported in part by grant DP0345001 from the Australian 
Research Council. This research has made use of the NASA/IPAC Extragalactic Database 
(NED) which is operated by the Jet Propulsion Laboratory, California Institute of Technol- 
ogy, under contract with the National Aeronautics and Space Administration. This research 
has made use of NASA's Astrophysics Data System. The Parkes telescope is part of the 
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Table 1: Resolution and Sensitivity 



9^02 


AS 


AT B 


AN HI 


S'(pred) 


AMpred) 


arscsec 2 


mJy/beam 


K 


10 18 cm- 2 


mJy 


10 6 M 


112x16 


0.57 


0.27 


8.2 


1.4 


4 


169x42 


0.61 


0.07 


2.2 


5.7 


15 


235x45 


0.71 


0.06 


1.8 


8.6 


23 



Notes: Columns give (1) half-power beam widths of synthesized beam, (2) rms noise level 
achieved in flux density in a single spectral channel, (3) rms sensitivity in brightness temper- 
ature, (4) rms sensitivity in HI columns density per 16.6 km s _1 velocity resolution, (5) the 
expected flux density if the beam were filled with Nhi = 10 19 ' 3 cm~ 2 emission in 16.6 km s _1 
, and (6) the HI mass contained in the beam, if filled with Nhi = 10 19,3 cm~ 2 . 
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Table 2: Nearby Galaxies 



Inn v"i t 

KjdbidbXy 


rvei. velocity 


rioj. oepaiation 




km s 1 


kpc 
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65 


V326 
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V376 
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V297 


104 


86 


V357 


1113 


93 


V397 


576 
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V294 




112 


W1216+064 


68 
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V327 


919 


123 


V278 


682 


129 


S1220+0638 


-399 


134 


V287 


179 


149 


NGC4241 


-1162 


155 


V332 


373 


156 


V384 




I6l 


V331 


-1238 


162 


S1219+0616 


-370 


172 


V436 




173 


V415 


665 


173 


V296 




174 


V349 




179 


S1218+0615 


595 


189 


V346 


-15 


190 



Notes: Galaxy names use V to indicate the VCC, S to signify an abbreviated SDSS entry and 
W for WSRT. The galaxy velocites are relative to the sub-DLA at 1895 km s _l . Projected 
separation of the sight line assumes a distance of 26.6 Mpc. 
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Fig. 1. — Image of the PG1216+069 field in HI 21cm emission over the heliocentric velocity 
range 900 to 2700 km s" 1 with synthesized beam 169" x42". Contours of beam-averaged, 
integrated HI brightness overlay the grayscale DSS optical image. Contours at 150, 300, 
450, 600, 750 mJy km s -1 per beam. PG1216+069 falls at the center. Labeled crosses 
indicate the galaxy centers for galaxies with redshifts between 895 and 2895 km s -1 , where 
the first part of the label is a galaxy ID (see Table 2) and the ±XXX indicates velocity 
relative to V hel = 1895 km s -1 . WHI indicates WSRT-J121918+0624 with velocity +68 
relative to the sub-DLA redshift. Triangles mark VCC galaxies having redshifts differing 
by more than 1000 km s" 1 from Vhel = 1895, while circles mark VCC galaxies of unknown 
redshift. X indicates NGC4223 (V he i = 2235), far outside the primary beam. Galaxies with 
V he i > 3000km s" 1 are omitted. 
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Fig. 2. — HI spectra for selected objects in the PG1216+069 field. Vertical dashed line 
indicates the absorption line redshift. Spectra for WSRT-J121921+0638 and VCC339 have 
been smoothed with a gaussian kernel with FWHM of 2 channels (16.6 km s -1 ); other spectra 
are unsmoothed. 



